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ABSTRACT 
A  Markov  Model  is  developed  to  study  the  distribution  of  faculty  in 
a  university.    The  purpose  of  this  paper  is  to  determine  the  underlying 
factors  that  control  the  distribution  of  faculty  members  from  one  time 
period  to  the  next.    Maintainable  and  realizeable  distributions  are  defined 
The  effects  of  appointment  policies  and  promotion  and  retirement  policies 
in  controlling  distributions  are  discussed.    The  model  is  tested  by 
comparing  the  predicted  output  with  real  data. 
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I.     INTRODUCTION 
In  any  graded  manpower  system,  one  of  the  principal  concerns  of 
top  level  management  is  the  control  of  the  distribution  of  personnel  among 
the  various  ranks  to  enhance  the  operation  of  the  organization.    The 
means  available  to  management  for  controlling  this  hierarchy  are  as 
follows: 

(a)  The  number  of  personnel  recruited  into  each  grade  at  each 
point  in  time . 

(b)  The  number  of  persons  promoted  (or  demoted)  from  one 
grade  to  the  next  within  the  organization. 

(c)  The  number  of  personnel  dismissed  from  the  organization 
or  induced  to  leave  the  system . 

In  most  organizations,  distribution  control  by  means  of  demotions  and 
dismissals  are  last  resort  measures  that  are  employed  only  under  extreme 
circumstances.    Control  by  (a)  alone  has  many  attractions  since  it  is 
least  disruptive  to  those  persons  already  in  the  system.     However,  con- 
trol of  unwanted  growth  at  the  top  is  usually  achieved  by  means  of 
reduced  promotion  rates  and  induced  retirement. 

The  purpose  of  this  paper  is  twofold:    first  to  design  a  mathematical 
model  that  quantitatively  describes  the  movement  of  university  faculty 
between  ranks  and  to  study  the  effects  that  promotion,  appointment  and 
retirement  policies  have  on  future  rank  structures  of  faculty;  and  second 
to  test  this  model  in  a  case  study  by  comparing  predicted  results  with 

actual  data.     It  has  been  recognized  for  some  time  that    promotion 
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practices  and  retirement  policies  have  a  profound  effect  on  future  faculty 
distributions.    Few  people  would  agree  on  an  ideal  distribution  of  faculty 
within  a  college  or  department,  although  most  persons  would  agree  that 
certain  distributions  are  undesireable .    For  instance,  a  distribution  of 
faculty  with  all  members  in  full  professor  rank  could  be  argued  as  unde- 
sirable both  from  an  economic  standpoint  and  because  such  a  distribution 
of  faculty  would  greatly  reduce  the  infusion  of  new  people  and  ideas  into 
the  institution.     No  attempt  will  be  made  in  this  paper  to  prescribe  a 
desired  distribution  of  faculty.    Rather,  this  paper  will  be  concerned 
with  a  quantitative  analysis  of  the  types  of  faculty  distributions  that  can 
be  realized,  given  certain  promotion  and  retirement  policies,  and  how 
these  distributions  will  be  effected  by  alterations  to  these  policies. 
The  following  assumptions  are  made  throughout  this  paper: 

(a)  The  organization  being  studied  (college,  department  of  a 
college,  or  university)  has  a  faculty  whose  total  number 
is  constant  from  year  to  year. 

(b)  The  movement  of  faculty  between  ranks  from  one  year  to 
the  next  is  a  Markov  type  process  in  that  it  is  possible  to 
deduce  the  future  distributions  from  knowledge  of  the  pre- 
sent rank  distribution.    Although  this  assumption  may  seem 
severe  and  in  some  circumstances  unrealistic,  the  results 
of  the  prediction  model  show  that  it  is  a  reasonable  one 

to  make  when  trying  to  predict  over  periods  of  5  years. 

This  paper  begins  with  a  general  description  of  the  model,  followed 

by  a  detailed  explanation  of  its  uses  in  studying  the  effects  of  policy 
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decisions.    This  is  followed  by  a  case  study  of  the  College  of  Engineering 
at  the  University  of  California,  Berkeley ,  where  actual  data  is  used  in 
estimating  the  Promotion  Matrix. 
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II,     THE  MARKOV  MODEL 

Consider  a  discrete  time  Markov  model  consisting  of  m  states  where 

state  j  is  the  rank  of  a  faculty  member  of  a  university ,  j  =  l,2...,m. 

These  m  states  can  be  divided  into  two  categories  namely:    n  transient 

states  and  m-n  absorbing  states.    The  sixteen  states  in  the  case  study 

are  defined  as  follows: 

States  Description 

1-4  Assistant  Professor  steps  1-4 

5-7  Associate  Professor  steps  1-3 

8-12  Full  Professor  steps  1-5 

13  Full  Professor,  overscale 

14  Retired  Professor 

15  Resigned  Professor 

16  Deceased  Professor        , 

where  states  1-13  are  the  n  =  13  transient  states  and  states  14  -  16 
are  the  m-n  =  3  absorbing  states. 

Let  P.,   be  the  fraction  of  faculty  in  state  j  on  1  July  of  a  particular 
year  that  were  in  state  i  as  of  1  July  of  the  previous  year.    The  date 
1  July  is  taken  for  administrative  bookkeeping  reasons.    Since  demotions 
are  extremely  rare  in  a  university  system,  we  take  P..    =  0  for  j  <  i.    The 
state  corresponding  to  retired,  resigned  and  deceased  faculty  are  con- 
sidered to  be  absorbing  states;  thus,  P..   =   1  for  i  =  14,   15,  and  16. 
Another  constraint  on  the  P . . '  s,  is  that  0  <  P. .   <  1    for  i  =  1 , . . . ,  13  which 
simply  requires  that  a  positive  fraction  of  the  faculty  in  any  given  state 
will  remain  in  that  state  for  the  succeeding  time  period,  and  that  no 
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active  faculty  member  can  remain  in  a  given  state  indefinitely.     Because 
the  state  space  describes  all  possible  states  of  nature  for  faculty,  this 
implies  that: 
m 


I     P..     :  1    ,      i  -  1,2 m 

i=i 


ij 


The  fractions,  P  . ,  are  conveniently  set  out  in  a  transition  matrix,  Q, 
as  follows: 


Pll         P13        Pi  3 

i 
I 

0  P22        P23 

0  0  P33 


Q     = 


Pi    i+i      Pi   n+s      Pi   m 


*2  n+i      ^2,n+2      ^2,m 
P3,n+i      P3,n+2      P3,m 


1 

0 

0 

0 

1 

o    ! 

0 

0 

i 

The  rows  of  the  Q  matrix  correspond  to  the  present  state  and  the  columns 

correspond  to  the  state  of  the  system  one  year  from  now. 

Consider  a  new  matrix  P  formed  from  Q  by  deleting  the  last  m-n  rows 

and  columns.    This  new  matrix,  called  the  Promotion  Matrix,   specifies 

the  transition  processes  of  active  faculty  members.     P  is  a  square  (  n  x 

n)  matrix  with: 

n 

I     P..    si,       i  =  l,2,.,.,n 

j=i  : 

and  at  least  one  row  of  the  P  matrix  sums  to  less  than  one. 
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Let  X*(t)  be  an  n-dimensional  row  vector,  X*(t)  =  (xf  (t) . .  .x*(t) )     , 
where  x*  (t)  represents  the  number  of  faculty  in  state  i  at  time  t.    Then 
X*(t)  is  a  vector  representation  of  the  distribution  of  faculty.    Further, 
let  the  appointments  in  year  t  be  represented  by  an  n-dimensional  row 
vector  (F*(t)  =(ff  (t)    . . .  f*(t)  )     ,  where  f *  (t)  is  the  number  of  faculty 
appointed  to  position  i  in  year  t  and  we  assume  that  there  is  no  con- 
straint on  the  availability  of  persons  to  fill  open    positions. 

The  distribution  of  faculty  at  time  t  +  1  is  determined  from  the  faculty 
distribution  at  time  t,  the  promotion  matrix,  and  the  appointment  policy 
at  time  t  +  1.    In  vector  matrix  notation  this  becomes: 

X*(t  +  1)  =X*(t)'P +  F*(t  +  1).  (1) 

Since  it  is  not  possible  to  have  a  negative  number  of  faculty  in  any 
rank  or  to  make  negative  appointments  to  a  rank  (dismissing  of  faculty 
is  considered  as  a  resignation  in  the  transition  matrix),  it  is  necessary 
that  x*(t)  and  f*(t)  *  0  for  i  =  1 ,2 , .  .  .  ,n 

If  the  total  faculty  size  is  to  remain  constant  from  one  year  to  the 

next,  the  total  number  of  personnel  leaving  the  system  each  year  must  be 

replaced  by  new  faculty.    Since  the  number  leaving  the  system  in  any 

given  year  is  a  function  of  the  faculty  distribution  during  the  year,  the 

total  number  of  positions  available  for  appointment  in  year  t  +  1  is  given 

by: 

n 

£    f*(t+  1)   =  X*(t).(I-P)a     ,  (2) 

i=l 

where  1    is  an  n-dimensional  column  vector  of  l's,  and  I  is  an  (n  x  n) 

identity  matrix. 

14 


III.    EFFECT  OF  APPOINTMENT  POLICIES 
Since  appointment  policies  are  the  most  desired  means  of  controlling 
faculty  distributions,  we  investigate  the  extent  to  which  distributions  can 
be  controlled  by  this  means .    Recall  that  the  total  number  of  faculty 
remains  constant  from  one  time  period  to  the  next.     Under  this  condition 
we  can  see  from  equation  2  that  the  number  of  positions  available  for  new 
appointments  in  a  given  year  is  a  function  of  the  previous  year's  faculty 
distribution.    In  order  to  better  understand  the  extent  to  which  distribu- 
tions are  restricted  by  the  transition  process,  three  subsets  of  the  set  of 
possible  distributions  are  defined  and  discussed,  namely: 

A)  The  Attainable  set  -  those  distributions  of  faculty  which  can 
be  attained  from  at  least  one  previous  distribution. 

B)  The  Maintainable  set  -  those  faculty  distributions  which  can 
be  maintained  from  one  time  period  to  the  next. 

C)  The  Containment  set  -  those  distributions  of  faculty  which 
can  be  attained  infinitely  often,  but  not  necessarily  in  consecutive 
periods.    This  is  followed  by  a  discussion  of  appointment  policies  and  how 
they  can  be  used  to  achieve  desired  distributions. 

In  studying  the  distribution  of  faculty  by  rank,  it  is  more  convenient 

to  consider  the  fraction  of  the  total  number  of  faculty  in  a  particular  rank 

rather  than  the  absolute  number  of  faculty  in  that  rank.    For  this  reason  we 

define  a  new  vector  X(t)  =  (xx(t), .  .  .xn(t)  ),  where 

x*(t) 
x,  (t)  = 


r  x*(t) 
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Note  that  in  the  remainder  of  this  paper,  if  no  limits  are  given  on  the  sum- 
mation sign  it  is  assumed  that  the  summation  is  over  all  transient  states 
and  the  dummy  variable  of  summation  is  omitted  when  no  ambiguity  exists  . 
Define  the  set, 

Y  =  { X(t)    |  X(t)  •  ,1  =   1 ,  xA  (t)  ^  0  for  all  i  }   . 
The  set  Y  then  contains  all  faculty  distributions. 

A.         THE  ATTAINABLE  SET 

The  set  of  all  distributions  is  the  fundamental  simplex  in  n-space 
(see  Figure  1).    The  question  naturally  arises,  do  there  exist  distributions 
to  which  we  can  return  in  some  future  year  but  which  might  not  be  main- 
tained from  year  to  year?    If  the  time  period  for  repetition  is  sufficiently 
small  and  the  intervening  distributions  sufficiently  close  to  the  desired 
distribution,  the  set  of  such  distributions  would  be  of  great  practical 
interest. 

A  necessary  condition  for  repeating  a  distribution  is  that  the  distri- 
bution must  be  reachable  from  at  least  one  distribution  in  Y.    The  set  of 
all  such  reachable  distributions  we  call  the  Attainable  set  A.    Since  it  is 
only  possible  to  make  non-negative  appointments,  F*(t)^  0  for  all  t,  and 
following  equation  1  we  define  A  as  follows: 

A={x|x^Z-P;X,ZeY} 
Although  this  definition  characterizes  the  Attainable  set,  it  does  not 
explicitly  describe  the  boundaries  of  this  set.    We  now  seek  a  character- 
ization of  A  which  describes  its  boundaries. 
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As  mentioned  above,  the  total  number  of  new  appointments  available 
in  any  year  is  a  function  of  the  distribution  of  faculty  the  previous  year. 
Define  a  column  vector,  W  =  (W1,  .  .  .  ,Wn)'  ,     as  follows: 

m 
W  =1        Pik      ,      i  =  l,2,...,n. 

k=n+l 

Then  W.  is  the  proportion  of  faculty  in  rank  i  which  annually  leave  the 

institution. 

Consider  the  extreme  points  of  the  set  of  all  distributions,  Y, 

represented  by  the  vectors     {  e.    }  ,  the  n-dimensional  unit  vectors,     Since 

the  extreme  points  are  contained  in  Y,  any  distributions  reachable  from 

the  extreme  points  must  be  elements  of  A.    Define  a  vector  F(t)  =  (^(t), 

.  .  .  ,fn(t)  )  where  f.  (t)  is  the  fraction  of  new  appointments  to  rani:  i  at  time 

t.     Then  equation  1  can  be  written  as  follows: 

X(t  +  1)  =  X(t)*P  +  X(t)-W-F(t  4   1)     .  (3) 

Let  X(0)  =  e.  .    Then  the  distributions  that  are  attainable  at  time  t  -  1  are 

given  by: 

X(l)    =    e    -P      +      e  -W-F(l)      ,  (4) 

i  i 

=    pW       +     Wt'F(l)      ,  (3) 

where  P^1'   is  the  ith  row  of  P  and  F(l)  is  any  appointment  vector. 

Consider  any  other  starting  point  Z  =  (Z^,  . .  .  ,  Zn) .    The  coordinates 

of  Z  can  be  written    LZ^e.,     Substituting  this  representation  of  Z  into 

equation  4  for  e^  we  get: 

X(l)      =      iZ.-e.-P      +     I  Z.-e.  'W-F(l)       , 
ii  ii 
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=      EZ.-P^      +      SZ.-W.  -F(l)       , 

=      EZ.  {   P(l)    +     W..F(1)    }      .  (6) 

Since  the  Z^s  sum  to  one,  any  attainable  distribution  given  by  equation 
6  is  a  convex  combination  of  the  distributions  described  by  equation  5. 
Thus  the  Attainable  set,  A,  is  the  convex  hull  of  the  points  with  co- 
ordinates: 

p(l)      +     W^e.    (j=l,2,...,n;  i=l,2,...,n)      (7) 
At  this  point  we  consider  an  example  to  illustrate  the  Attainable  set 
In  order  to  illustrate  the  problem  diagramatically,  we  consider  the 
following  reduced  state  space: 

State  1  =    Assistant  Professors 

2  =  Associate  Professors 

3  =  Full  Professors 

4  =  Retired,  resigned  or  deceased  Professors. 
A  typical  transition  matrix  may  be  as  follows: 

12  3  4 


Q 


.5 

.4 

0 

.1 

0 

.6 

.3 

.1 

0 

0 

.7 

.3 

0 

0 

0 

1. 

The  promotion  matrix  would  then  be: 


P     = 


.5 

.4 

.0 

.0 

.6 

.3 

.0 

.0 

.7 
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ASSOC, 
PROF 


(0,  5,.5) 


(0,0,1) 

FULL 
PROF 


(.3,0,7) 


(.5,  .5,0) 

(.6,-4,0) 


FIGURE  I 


THE  ATTAINABLE  SET 
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Equation  7  gives  the  following  points  of  the  Attainable  set: 

(.5  +  .1,  .4,  0)*,  (.5,  .4  +  .1,  0)*,  (.5,  .4,  0  +  .1), 
(0  +  .1,  .6,  .3),  (0,  .6  +  .1,  .3)*,  (0,  .6,  .3  +  .1), 
(0  +  .3,  0,   .7)*,        (0,  0  +  .3,   .7),         (0,  0,   .7  +  .3)*     . 

The  extreme  points  of  A  are  marked  with  an  asterisk.    The  boundaries  of 

this  set  are  illustrated  in  Figure  1 . 

B.         THE  MAINTAINABLE  SET 

The  set  of  all  faculty  distributions  with  the  property  that  the  distri- 
bution can  be  duplicated  in  consecutive  time  periods  is  called  the  Main- 
tainable set  and  is  denoted  M.    In  the  literature,  this  is  also  called  the 
set  of  steady  state  distributions  (see  Bartholomew  (1967)  ). 

Let  F(t)  =  F  for  all  t;  that  is  we  appoint  the  same  distribution  of 
faculty  each  year.    From  equation  3  the  steady  state  distribution  of 
faculty,  X,  is  given  by:    X  =  X*P  +  X»W'F  .    Let  F*(X)  =  X-W-F  .    ThenF* 
gives  the  number  of    faculty  which  must  be  appointed  to  each  rank  in 
steady  state  to  maintain  the  distribution  X.    It  follows  that  F*  =  X*  (I-P). 
Since  F*  must  be  non-negative,  we  now  define  M  as  follows: 

M  =  {  X  |  X-(I-P)  £  0,  X  c  Y}   . 
Since  at  least  one  row  of  the  P  matrix  sums  to  less  than  one,  the  matrix 
(I-P)  has  an  inverse  and  all  terms  of  this  inverse  are  greater  than  or  equal 
to  0  (see  Kemeny  and  Snell,     (1960)  ),     Define  N  =  (I-P)-1.    Then 

X  =  F*'N     .  (9) 
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In  what  follows  it  is  more  convenient  to  consider  the  fraction  of  new 
appointments  to  each  particular  rank  rather  than  the  total  number.    For 
this  reason  we  define  a  new  matrix,  B  =(B..),  formed  from  the  N  matrix  by 
normalizing  the  row  sums;  hence  let 


N.. 
ij 


B 


ij 


N 


k 


ik 


(9) 


The  necessary  condition  for  X  to  be  a  member  of  the  Maintainable 
set  as  stated  in  equation  8  can  now  be  written: 

X  =  F-B    ,  where  (10) 

f  .      nf*'EM, 

i         i      .       ik 
k 

Since  we  require  that  EX.    =    1 ,  it  follows  from  the  definition  of  B 

that  Lf.  =  1,  the  desired  result.    Since  f.  ^  0   and  If.  =  1,  equation  10 
l  11 

expresses  X  as  a  convex  linear  combination  of  the  rows  of  B.    Thus  the 
rows  of  B  represent  the  extreme  points  of  the  Maintainable  set  which  is  a 
closed  convex  polyhedron  in  n-space. 

Referring  to  the  example  of  Section  III  A,  we  compute: 


N=      (I-P)-1 


N. 


B  = 


ij 


L  N 


ik 


2.0 

2.0 

2.0 

0 

2.5 

2.5 

0 

0 

3.3 

.33 

.33 

.33 

0 

.50 

.50 

0 

0 

1.0 

The  Maintainable  set  is  the  indicated  region  of  Figure  2.    The  extreme 
points  of  this  Maintainable  set  are  described  by  the  three  vectors  defined 
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ASSOC. 
PROF 


(0.5..5) 


1:3  NON-TENURE /TENURE 
FACULTY  RATIO 


rMAINTAINABLI 

SET^ 


,(.33,. 33,  .33) 

REDUCED  SET  DUE  TO 
_.  N£)NrTEW«£-C.QN5IRAINT. 


(0,0,1) 
FULL 
PROF. 


FIGURE    2. 
THE  MAINTAINABLE  SET 
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by  the  rows  of  B,  namely  (.333,    .333,   .333),     (0,    .50,    .50)  and  (0,  0,   1). 
As  can  be  seen  from  Figure  2,  the  Maintainable  set  is  only  a  small  portion 
of  the  set  of  all  distributions  with  a  tendency  toward  "top  heavy"  faculty 
distributions.    If  a  constraint  is  added  specifying  that  the  ratio  of  non- 
tenure faculty  be  no  less  than  1:3,  the  constrained  Maintainable  set  is 
considerably  reduced  as  illustrated  by  the  cross-hatched  region  of  Figure 
2.     It  is  easy  to  see  that  many  non-tenure/tenure  ratios  are  not 
maintainable. 

C.         THE  CONTAINMENT  SET 

The  primary  reason  for  defining  the  Attainable  set  was  to  restrict  the 
set  of  all  distributions  to  those  which  can  be  reached  at  least  once  from 
any  current  distribution.    We  now  wish  to  find  the  set  of  distributions 
which  can  be  reached  infinitely  often,  though  not  necessarily  in  successive 
time  periods.    We  shall  call  this  new  set  the  Containment  set  and  desig- 
nate it  by  the  letter  C  .    In  Section  III    A  it  was  shown  in  equation  6  that 
the  distributions  reachable  from  any  given  starting  distribution  are  a  sub- 
set of  the  Attainable  set.    Also,  the  Maintainable  set  is  a  subset  of  the 
Containment  set  since  any  distribution  in  the  Maintainable  set  can  be 
achieved  infinitely  often  with  a  constant  appointment  policy.     JJence  we 
have  M  a  C   «=A.     Mathematically  we  define  C  as  follows: 
C   =    [X    I  X  e  Y,  X  €  Rn(X)  for  some  n  }  ,  where 
Rn(X)    =    {X  |  X,X  cY,  X  is  reachable  from  X  in  n  steps   }    , 

=    {  X   |  X  ,X  €  Y ,  X  =  X«  Pn      +     I  f  *  (i) '  P0"1  } 

i 

for  some  sequence  of  appointment  vectors  (f*(i)}  . 
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Although  we  have  defined  C,  the  definition  provides  little  insight 
in  describing  the  boundaries  of  C.    Unlike  the  boundaries  of  the  Maintain- 
able and  Attainable  sets  we  are  unable  to  provide  a  precise  description  of 
the  boundaries  of  the  Containment  set.    We  conjecture  that  these  bound- 
aries form  an  infinite  sided  polyhedron.    Some  calculations  in  the  3-dimen- 
sional  case  indicate  this.    We  proceed  at  this  point  to  obtain  approxima- 
tions of  C  in  the  3-dimensional  case  to  gain  some  idea  of  its  size  and  its 
relation  to  the  Maintainable  set. 

It  is  apparent  that  any  distribution  lying  outside  the  Maintainable 
set  that  is  reachable  from  a  point  within  the  Maintainable  set  in  n  steps 
must  be  an  element  of  C  since  return  to  the  starting  distribution  is 
guaranteed  by  definition  of  M.    Using  this  information,  a  computer  pro- 
gram was  developed  to  delineate  those  distributions  reachable  from  points 
of  M  using  appointment  policies  designed  to  keep  the  resulting  distribu- 
tions out  of  the  Maintainable  set.    In  the  example  of  Section  III    A,  the 
distributions  resulting  from  appointing  only  Assistant  Professors  to  an 
initial  distribution  of  only  Full  Professors  were  computed.    This  proceedure 
was  continued  until  the  number  of  Associate  Professors  started  to  decrease 
in  number.    At  that  point,  the  appointment  policy  was  changed  to  a  policy 
of  appointing  only  Associate  Professors.    It  is  conjectured  that  in  the  3- 
dimensional  case  this  approach  gives  the  boundaries  of  C  as  displayed  in 
Figure  3.    It  is  not  clear  how  this  construction  would  be  extended  to  higher 
dimensions. 
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Another  approach  to  approximating  the  boundary  of  C  in  the  3-dimen- 
sional  case  is  by  means  of  constraint  equations  on  the  various  ranks. 
Starting  with  any  distribution  of  faculty  in  the  Attainable  set  at  time  t  we 
maximize  the  number  of  Assistant  Professors  in  time  t+1  subject  to  the 
constraint  that  the  number  of  Assistant  Professors  does  not  decrease  in 
years  t  to  t+1.    We  can  see  from  Figure  3  that  such  a    proceedure  provides  a 
an  upper  bound  on  the  Assistant  Professor  coordinate  of  the  Containment 
set.    A  similar  proceedure  was  used  to  bound  the  Associate  Professor 
coordinate  of  C .    To  find  a  lower  bound  for  the  Full  Professor  co-ordinate, 
X3,  we  maximized  Xx   +  X2  subject  to  the  constraint  that  this  sum  does  not 
decrease  from  time  t  to  t+1. 

Using  the  example  of  Section  III    Af  the  following  calculations  are 
displayed  in  Figure  3. 

1.  Maximize  X^t+1)  =  X^tJ'Pu  +  X(t)'W-e1 

=  X1(t)-(P11   +  P14)   +  X^tJ'P^  +  X3(t)-P34 

=  .GX^t)   +   .lX2(t)    +    .3X3  (t) 
Xx(t+1)    ^  Xx(t)   ^X^t)    £  .25X2(t)    +    .75X3(t) 
X(t)€Y  =*  Xa(t)  =  1   -  Xa(t)    -  X2(t) 

;.       Xx(t)      £     3/7   -   2/7X2(t) 
X(t)  e  A  =»    Maximum  X^t+1)   =    .393 

2.  Maximize  X2(t+1)   =  X^-P^  +  X2(t)  •  P^  +  X(t)-W.e2 

=  Xa  (t) .  (Pia  +  P14 )  +  X2(t)  •  (Pss  +  Ps, )  +  X3(t)  •  P34 
=    .SX^t)  +  .7X2(t)  +  .3X3(t) 
X2(t+1)  ^  X2(t)    =>X2(t)  £  5/3Xx(t)   +  X3(t) 


26 


X(t)  €  Y     «    X^t)   =   1    -  X2(t)    -  X3(t) 

.'•       X2(t)    <    5/8   -   2/8X3(1:) 
X(t)    c  A     =»    Maximum  X8(t+1)   =    .591 
3  .    Maximize  Xx(t+1)  +  X2(t+1)  =  Xx(t)  +  Xs(t)  *  (P22  +  PaJ  +  X3(t)  •  P34 

=  X1(t)  +  .7X2(t)  +  .3X3(1:) 
X^t+1)  +Xs(t+1)  ^Xx(t)  +X2(t)  =>X3(t)    *Xa(t) 

X(t)  e  A  and  Maximum  X^t+1)  =  .393  =>  Maximum  Xx(t+1)  +  X2(t+1) 

=  .607 

D.        APPOINTMENT  POLICIES 

Up  to  now  we  have  discussed  the  set  of  distributions  attainable 
under  various  appointment  policies.     In  this  section  we  discuss  the  types 
of  appointment  policies  necessary  to  achieve  certain  faculty  distributions. 
Since  distributions  in  the  Maintainable  set  can  be  maintained  by  means  of 
a  constant  appointment  policy  as  shown  in  equation  10,  it  will  be  assumed 
that  the  desired  distribution  lies  outside  the  Maintainable  set  but  within 
the  Containment  set.    Since  distributions  in  this  region  are  not  maintain- 
able from  time  period  to  time  period,  the  sequence  of  distributions 
realized  during  a  repeating  cycle  contains  elements  of  both  the  Contain- 
ment set  and  the  Maintainable  set.    We  limit  our  study  of  appointment 
policies  to  considering  their  effect  on  distributions  in  the  Maintainable 
set  since  these  effects  are  easier  to  analyze.    Let  X  be  the  steady  state 
distribution  of  faculty  such  that  X.    =  fraction  of    total  faculty  in  rank  i. 
Recall  from  equation  10  that  X  =  F«B. 
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Due  to  the  upper  triangular  feature  of  the  promotion  matrix  and  the 
resulting  similar  configuration  of  the  N  and  B  matrices,  the  effect  of 
making  all  new  appointments  at  one  level  (i.e.  making  F  a  unit  vector) 
will  realize  a  distribution  of  faculty  as  described  by  the  row  of  the  B 
matrix  corresponding  to  the  level  of  appointments.    This  means  that  a 
certain  faculty  rank  cannot  be  maintained  at  a  zero  level ,  or  the  probabil- 
ity of  promotion  is  zero.    Mathematical^ ,  this  becomes: 

xi  a  7^-|Ixi-pij|        'andso 

X.   =  0   =»    either  X.   =  0        i  <  j 
or  P . .   =  0  i  <  j   . 

In  order  to  simplify  analysis  of  the  model  we  rule  out  the  possibility  of 
double  and  triple  promotions.    With  this  simplification,  the  above 
inequality  reduces  to: 

Xj    *     1-P  *  W'J       ' 

If  the  probability  of  remaining  in  a  particular  rank  plus  the  probabil- 
ity of  being  promoted  into  that  rank  sums  to  more  than  one,  then  the  number 
of  faculty  in  that  rank  must  be  greater  than  the  number  in  the  preceding 
rank.     Mathematically,  if 

Pi!    j   +   P.  ,,    >  1    =»   Xi    >    X,  .,     and  if 

pi-i,j  +  p)j  -  1=>xi  a  xj-i  • 

Since  a  rank  can  only  be  maintained  by  appointing  at  that  rank  or  at 
lower  ranks,  it  is  necessary  to  make  some  portion  of  the  total  appoint- 
ments at  the  lowest  rank  in  order  that  all  ranks  be  represented.     This  point 

is  obvious  if  one  stops  to  consider  the  implications  of  the  promotion  matrix. 
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Consider  the  sixteen  states  defined  on  page  4.    For  the  largest 

ratio  of  non-tenure  faculty  to  total  faculty  (where  non-tenure  faculty 

includes  Assistant  Professor,  steps  1-4),  appointments  should  be  made 

exclusively  at  that  rank  where  the  sum  of  the  first  four  elements  of  the  B 

matrix  is  maximum.    The  above  appointment  policy  statement  can  thus  be 

expressed  as  follows:      appoint  to  state  i,  where  i  is  chosen  so  that 

4  4 

I     B^      £    I    Bkj     for  all  k  =  l,2,...,n  . 

j=l  j=l 

For  the  case  where  only  single  promotion  advancements  are  permitted, 

appointments  should  be  made  exclusively  at  the  Assistant  Professor  Step 

1  level  since: 


IB.       ^    I     B,  .      i  <  k 
.    ,    U  .    ,    kj 

j  =  l  j  =  l 


Proof: 


13 
j  =  l 


£    B.       =1    since  B  is  normalized. 


13 


E    B. ,      =      1  -    £  B.       for  all  i  =  1,2,  .  .  .  ,n 
j=l  j=5 


For  single  promotion  only,  it  can  be  shown  that 

p.    . 
ri,i+i 

B 


i/i 


If  Pu   +   Pi/i+1      =      1      *     Bij    =   Bi+1/jl    and  if 
pii   +  pi,i.i     <      1      "     Bij    <   Bi+1J    . 
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Therefore,  13  13  and  hence 

j=5  j=5 


I   Bij     >     L    Bkj  i  <  k     . 


4  4 

E    B^     *     £ 
J-l  J-l 

Thus,  one  can  increase  the  non-tenure/tenure  ratio  by  appointing 
at  the  Assistant  Professor  Step  1  level  and  reduce  this  ratio  by  appointing 
at  the  Assistant  Professor  Step  4  level.     Such  effects  of  appointment 
policy  are,  however,  predicated  on  the  fact  that  the  transition  matrix  will 
remain  essentially  unchanged.    The  effect  of  changes  in  the  transition 
matrix  will  be  discussed  in  the  next  chapter. 

The  effect  of  a  mixed  appointment  policy  is  easily  determined  once 
the  B  matrix  is  known.    A  policy  of  appointing  20%  of  the  new  appoint- 
ments to  tenure  positions  could  be  evaluated  by  taking  80%  of  row  1  of  the 
B  matrix  (corresponding  to  appointment  at  the  Assistant  Professor  Step  1 
level)  and  20%  of  row  5  of  the  B  matrix  (corresponding  to  appointments  at 
the  Associate  Professor  Step  1  level).    Such  an  appointment  schedule  will 
realize  the  greatest  achievable  spread  among  the  various  ranks,  given  the 
constrained  appointment  policy.    Other  appointment  policies  can  be 
evaluated  in  a  similar  way. 
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IV,     EFFECT  OF  PROMOTION  POLICIES 
Although  the  preferred  method  of  regulating  faculty  distributions  is 
by  means  of  an  appointment  policy  ,  it  is  evident  from  the  preceding 
chapter  that  other  means  must  be  employed  if  certain  distributions  are 
to  be  realized.     One  such  method  involves  changing  the  promotion  and/ 
or  retirement  policies  of  the  institution.    For  example,  if  a  general  policy 
of  requiring  faculty  members  to  remain  a  minimum  of  two  years  in  each 
grade  were  altered  to  a  three  year  minimum ,  the  change  would  have  a 
noticeable  effect  on  the  transition  matrix  and  subsequently  on  the  bound- 
aries of  the  Maintainable,  Attainable  and  Containment  sets.    Although 
the  exact  effect  that  policy  changes  have  on  transition  frequencies  is 
not  known,  these  relationships  can  be  determined  with  reasonable 
accuracy  by  analysing  past  transition  matrices  both  before  and  after  such 
changes  are  implemented. 

In  this  section,  we  examine  how  changes  to  the  promotion  matrix 
affect  maintainable  faculty  distributions  by  first  studying  individual 
transition  frequency  changes  and  then  looking  at  changes  to  all  the  tran- 
sition frequencies  associated  with  a  given  state.    In  order  to  simplify  the 
analysis  and  to  focus  attention  on  major  effects,  we  consider  only  single 
promotions  as  a  means  of  advancement.    Although  jump  promotions  are 
possible,  their  relative  occurrence  is  small  and  do  not  appreciably  alter 
the  results  of  the  analysis. 

With  the  above  restriction,  the  elements  of  the  N  matrix  (=  (I-P)-1 ) 

can  be  expressed  as  follows: 
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N        =  -^-       n        Pk,k+1  f or  j   ^  i    , 

J  Fjj       k=i      p* 

k,k 

=  0  for  j  <  i    , 


where      P*      =  1   -  P-      ,      and 

P*      =  P. . 
This  can  be  rewritten  in  recursive  form: 


p  * 
Nu  =    N,    K  '     HJ forj  ^i    , 

j  .j 


p*. 
JJ 


for  j  =  i,    and 


=    0  otherwise. 

In  this  form  it  is  more  apparent  what  effect  a  change  in  promotion  policy 
has  on  the  distribution  of  faculty. 

If  the  probability  of  remaining  in  state  j  ,  given  that  one  starts  in 
state  j,  is  increased  {?■■    new  >  P.M   old),  then  the  proportion  of  faculty 
in  states  j  through  n  will  increase  and  the  ratio  of  faculty  in  states  1 
through  j-1  will  decrease,  given  that  some  appointments  are  made  at  the 
j-1  level  or  below.    This  is  because  an  increase  in  P.:  j   results  in  a 

decrease  of  1  -P:  ■  or  an  increase  in  - — r —  .    Therefore,  in  the  N  matrix, 

J  J  i  -  r ,  • 

j) 

all  states  from  j  through  n  will  be  increased  proportional  to 


1-pjj 


while  states  1  through  j-1  will  remain  unchanged;  this  results  in  changed 
ratios  of  faculty  in  the  different  ranks.    The  opposite  result  holds  if   P... 
is  decreased. 
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If  the  probability  of  being  promoted  to  state  j+1,  given  one  starts 
in  state  j ,  is  increased,  the  ratio  of  faculty  in  states  j+1  through  n  will 
increase,  and  the  ratio  of  faculty  in  states  1  through  j  will  decrease, 

given  that  some  appointments  are  made  at  the  j  level  or  below.    This  is 

p 
because  an  increase  in  P..  results  in  an  increase  in  ij  .    Thus 

lj  1-P.. 

JJ 

in  the  N  matrix,  all  states  from  j+1  through  n  will  be  increased  by  an 

P 
amount  ij  ,    while  states  1  through  j  will  remain  unchanged.    The 

1-P.. 
JJ 

opposite  result  again  holds  for  a  decrease  in  P.. . 

When  both  the  probability  of  promotion  P.    .    ,  ,  and  the  probability  of 

remaining  in  state  j,P,    ,,  are  increased,  the  overall  effect  is  to  increase 

the  ratio  of  faculty  in  states  j+1  through  n,  given  that  appointments  are 

made  at  the  j  level  or  below.    The  effect  on  the  ratio  of  faculty  in  state  j 

is  dependent  upon  the  relative  increases  in  probabilities  and  upon  the 

appointment  policy.    If  appointments  are  limited  to  state  j  and  above,  the 

ratio  of  faculty  in  state  j  will  decrease.    If  the  increase  in  P..  is  much 

JJ 

greater  than  the  increase  in  P.    .      ,  and  appointments  are  made  at  the 

J  /  J    -*- 

state  j-1  level  or  below,  then  some  increase  in  the  ratio  of  faculty  for 
state  j  may  be  realized,  but  at  no  time  will  the  increase  be  as  large  for 
state  j  as  for  higher  states. 

When  the  probability  of  remaining  in  state  j  is  increased  while  the 
probability  of  being  promoted  from  state  j  is  decreased,  the  overall  effect 
on  faculty  distribution  is  to  increase  the  proportion  of  faculty  in  state  j 
and  decrease  the  proportion  of  faculty  in  states  j+1  through  n .    The  effect 
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on  states  1  through  j-1  is  dependent  upon  the  size  of  the  increase  or 
decrease  relative  to  the  initial  probabilities  and  the  subsequent 
appointment  policy. 
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V.    A  CASE  STUDY 
The  following  sections  describe  the  data,  transition  matrix,  and  a 
comparison  of  predicted  and  actual  distributions  resulting  from  a  case 
study  of  the  College  of  Engineering  at  the  University  of  California  at 
Berkeley. 

A.      COLLECTION  OF  DATA 

The  data  on  faculty  movements  within  the  College  of  Engineering  was 
gathered  by  examining  College  and  University  records  for  discrete  time 
points  in  the  past,  namely  at  12:01  A.M.  on  1  July  of  the  years  1960 
through  1968.    Faculty  in  a  particular  state  at  these  points  in  time  were 
considered  to  remain  in  that  state  throughout  the  one  year  time  period. 
The  data  was  collected  in  the  following  manner: 

1.  Historical  records  of  all  presently  employed  faculty  members  of 
the  College  of  Engineering  were  obtained  by  the  Dean's  Office. 
The  present  age  and  time  spent  in  each  position  was  recorded 
for  each  faculty  member,  and  accumulated  for  each  department 
in  the  college.    This  information  was  then  used  to  form  a  histor- 
ical record  of  the  faculty  distribution  of  the  various  departments, 
and  for  the  College  of  Engineering  as  a  whole,  from  19  60 
through  1968.    At  no  time  were  confidential  personnel  records 

or  biographical  data  made  available  to  us. 

2.  In  conjunction  with  the  above,  data  was  also  collected  at  Uni- 
versity Hall  to  check  the  lists  of  faculty  of  the  College  of 
Engineering  for  the  years  1960  through  19  68.    Starting  in  I9  60, 
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the  listing  of  faculty  members  by  department  was  recorded ,  and 
this  list  was  updated  through  19  68.    The  record  of  faculty 
members  currently  holding  positions  in  the  College  of  Engineer- 
ing was  compared  with  the  listing  from  the  Dean's  Office  to 
provide  an  independent  check  of  promotion  data.    In  addition,  a 
listing  was  made  of  faculty  members  who  were  employed  some- 
time during  the  period  1960  to  1967  but  who  are  not  presently 
employed,  in  order  to  provide  a  historical  record  of  the  entire 
College  during  this  period. 

3.  A  separate  listing  was  provided  by  the  Dean's  Office  recording 
the  total  numbers  and  distribution  of  faculty  from  1960  through 
1968,  as  ascertained  from  college  records.    Included  with  this 
list  was  a  breakdown  by  name  of  those  faculty  members  no 
longer  employed  by  the  University,  and  why  and  when  they 
terminated  employment.    This  information  was  used  as  a  check 
against  the  other  two  sources  of  information.     Only  in  a  few 
special  cases  was  the  data  from  these  sources  in  conflict; 
these  discrepancies  were  resolved  by  corresponding  directly 
with  the  Dean's  Office. 

4.  A  final  source  of  information  of  the  College  of  Engineering  was 
provided  by  an  independent  study  conducted  by  the  Chancelor's 
Office.     Using  information  gathered  from  University  files  as 
input  data,  a  computer  program  designed  by  personnel  of  the 
Chancellor's  Office  was  used  to  provide  a  complete  profile  of 
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the  College  of  Engineering  by  age  and  position  for  each  year 
from  1960  through  1968.    Although  this  information  was  gathered 
from  different  records  than  those  used  in  this  report,  and  a 
different  time  base  was  used,  the  overall  results  of  this  report 
agree  quite  closely  with  this  profile  listing. 
A  compilation  of  the  raw  data  is  presented  in  Tables  1  through  8. 

The  cumulative  data  for  the  period  1960  through  1968  is  presented  in 

Table  9. 

B .      THE  TRANSITION  MATRIX 

The  transition  matrix  for  the  College  of  Engineering  was  calculated 
in  the  following  manner: 

1 .  The  total  number  of  faculty  in  each  position  for  a  particular 
year  was  determined  by  summing  the  statistics  of  the  individual 
departments  within  the  College. 

2.  In  a  similar  manner,  the  total  number  of  faculty  in  each  position 
that  were  given  merit  increases,  promoted  or  accelerated  as 
well  as  those  who  resigned,  retired,  or  died  during  a  particular 
year,  were  recorded. 

3.  In  order  to  provide  a  larger  sample  base  for  estimating  para- 
meters, the  annual  statistics  for  the  period  1960  through  1968 
were  aggregated  to  form  totals  representing  a  single  period. 
This  procedure  is  justified  if  the  transition  process  is  in  fact 
Markovian,  and  if  no  major  policy  changes  occurred  during 


37 


this  time  period  that  could  effect  the  transition  process.    This 
author  was  able  to  discover  only  one  policy  change  during  the 
period  that  might  have  affected  the  transition  process:  the 
addition  of  new  steps  to  the  Assistant  Professor  and  Full  Pro- 
fessor ranks  in  the  period  1962-1963.    A  straightforward  aggre- 
gation of  the  data  from  1960-19  68  would  indicate  an  unrealistic 
double  promotion  rate  from  Assistant  Professor,  Step  3  to 
Associate  Professor,  Step  1,  since  the  Assistant  Professor, 
Step  4  state  did  not  exist  from  1960-1963.    We  are  not  concerned 
in  this  paper  with  the  actual  policies  used  by  the  College  of 
Engineering  during  the  transition  period.    Rather,  we  must 
devise  some  method  to  estimate  current  realistic  transition 
frequencies,  P33,  P34,  P^,  and  P45  from  past  data.    We  accom- 
plished this  by  considering  all  Assistant  Professors,  Step  3  to 
be  Assistant  Professors,  Step  4  in  the  period  1960-1963.    Thus, 
all  promotions  from  Assistant  Professor,  Step  3  to  Associate 
Professors,  Step  1  in  this  period  were  single  promotions.    Also, 
in  the  period  1960-1963,  since  Full  Professor,  Step  5,  did  not 
exist,  straightforward  aggregation  of  data  would  lead  to  un- 
realistic estimates  of  current  promotion  rates  from  Full  Profes- 
sor, Step  4,  to  the  overscale  level.    By  assuming  all  Full 
Professors,  Step  5,  were  overscale  in  1960-1963  we  were  able 
to  obtain  more  realistic  transition  frequencies  from  Full  Pro- 
fessor, Step  4  to  Full  Professor,  Step  5. 
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4.      The  elements  of  our  transition  matrices  were  derived  by  deter- 
mining the  observed  relative  frequency  of  transitions  between 
states.    Thus,  the  fraction  of  Assistant  Professors,  Step  1, 
that  become  Assistant  Professor,  Step  2,  in  the  next  year  is 
equal  to  the  observed  number  of  promotions  from  state  1  to  state 
2  divided  by  the  total  number  in  state  1  at  the  beginning  of  the 
year. 

The  transition  matrix  based  upon  the  data  from  the  years  1960  through 

1965  is  shown  in  Figure  4. 

C.     COMPARISON  OF  PREDICTED  AND  ACTUAL  DISTRIBUTIONS 

In  order  to  determine  how  well  the  model  predicts  actual  faculty 
distributions,  data  from  the  years  1960  through  1965  were  used  to  deter- 
mine transition  estimates  for  the  promotion  matrix.     Since  these  transition 
estimates  agree  quite  closely  with  the  estimates  based  upon  the  nine 
year  period  19  60  through  1968,  the  Markov  assumption  appears  to  be 
reasonably  well  justified.     The  distribution  of  faculty  in  1  July  1965  was 
used  to  calculate  returning  faculty.    New  appointments  were  added  to 
predict  the  distribution  of  faculty  as  of  1  July  19  66.    This  procedure  was 
repeated  for  subsequent  years,  using  the  previous  years  predicted  dis- 
tribution plus  the  actual  appointment  policy  of  the  current  year  to 
determine  the  current  distribution  of  faculty.    The  transition  matrix  based 
upon  the  data  from  the  years  1960  through  1965  is  shown  in  Figure  5. 
The  predicted  distribution  of  faculty  for  the  years  1966  -  1968,  based 

39 


upon  the  1965  distribution  of  faculty  and  the  ensuing  appointment 
schedules,  is  shown  in  Figure  6  along  with  the  actual  distributions  for 
these  years.    In  a  three-year  average,  the  estimate  of  Assistant  Profes 
sors  was  4.25%  below  the  actual,  the  estimate  of  Associate  Professors 
was  10%  above  the  actual,  the  estimate  of  Full  Professors  was  4.9% 
below  the  actual.    The  estimated  total  number  of  faculty  for  each  year 
was,  on  the  average,   1.7%  below  the  actual  total  number  of  faculty. 
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